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Abstract 

Much effort has been invested on effective-field-theoretical studies of the 
near-threshold NN^NNtt reactions and, in order to deal with the somewhat 
large three-momentum transfers involved, the momentum counting scheme (MCS) 
was proposed as an alternative to the usual Weinberg counting scheme. Given 
the fact that a quantitative explanation of the existing high-precision NN^ 
NNtt data requires a careful examination of higher chiral order contributions 
to the transition operator, we make a detailed numerical investigation of the 
convergence property of MCS for a pilot case of the pp — > pptr reaction. Our 
study indicates that MCS is superior to the Weinberg scheme in identifying 
dominant higher order contributions to the NN^NNir reactions. 



There exists a substantial accumulation of high-precision data on various observ- 
ables for the near-threshold NN^NNtt reactions, e.g. [I], and providing a coher- 
ent understanding of these experimental results has been a prominent theoretical 
challenge, see e.g., [21ElffllSllE10mi3]- Heavy-baryon chiral perturbation theory 
(HB;\PT), which is a low-energy effective field theory of QCD, is expected to offer a 
systematic framework for describing these reactions. In HB;\PT, the four-momentum 
Q characterizing a given physical process is assumed to be small compared to the 
chiral scale A x ~ Anf^ ~ 1 GeV, and contributions to the transition amplitude are 
classified according to the power (chiral order) of the expansion parameter e = (5/A x Q 
The coefficients of possible terms in the HB%PT Lagrangian, called low-energy con- 
stants (LECs), can in principle be linked to the matrix elements of QCD operators 
but in practice they are determined from experimental observables. Once all the LEC 
are fixed up to a specified chiral order, HB^PT allows us to make predictions for a 
wide range of hadronic and electroweak processes. 

Although, as mentioned, HB^PT presupposes the smallness of its expansion pa- 
rameter Q/A x , the pion production reactions involve somewhat large momentum 
transfers, p ~ y/m n mN, even at threshold. This implies that the application of 
HB%PT to the NN — > NNn reactions may involve some delicate aspects, but this 
also means that these processes may serve as a good test case for probing the ap- 
plicability (or the limit of applicability) of HB^PT. It is worth emphasizing that 
finding a valid EFT expansion scheme for the NN — >NNtt reactions is of general 
importance (going beyond the specific context of the NN^NNtt reactions) because, 
once such a scheme is found, we expect to be able to develop similar frameworks for 
other nuclear processes which involve rather large energy-momentum transfers and 
hence may not be quite amenable to the straightforward application of the ordinary 
HBxPT approach. To account for the rather large momentum transfers involved in 
the NN^-NNn reactions, Cohen, Friar, Miller and van Kolck [5] proposed to re- 
place the ordinary chiral counting scheme of Weinberg (W-scheme) with a new scheme 
called the momentum counting scheme (MCS). The expansion parameter in MCS is 
e = p/rriN ~ (m n /miy) 1 ^' 2 ~ 1/2.6, which is appreciably larger than the expansion 
parameter e ~ m^jm^ ~ 1/6.7 in the W-scheme. 

Although the formal aspects of MCS have been discussed rather extensively [5], 
it seems fair to say that its practical utility is yet to be established. An important 
point to be noted here is that, whereas in the W-scheme a given Feynman diagram 

1 The numerical value of the chiral expansion parameter is e ~ m T /A x ~ 1/7.1, while the "recoil- 
correction" expansion parameter is m^/rriN — 1/6.7. Since these parameters have roughly the same 
numerical value, the chiral and recoil-correction expansions are combined together in HB^PT. 
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corresponds to a definite power in e (it thus has a unique HB%PT order), this type of 
correspondence does not in general exist in MCS, because a given Feynman diagram 
can involve contributions that belong to different orders in e. In the following, the 
contribution corresponding to the lowest power of e for a given diagram is referred to 
as the "leading term" [8], and the remaining (higher order in e) contributions of 
the diagram as "sub- leading terms". It is to be emphasized that, in order to study 
the convergence property of MCS, we need to examine not only the behavior of the 
leading terms but also that of the sub-leading terms. Our first attempt at such a study 
was described in Ref. [10J [11], to be referred to as KSMK1. The results presented in 
KSMK1 indicate that the individual terms in the MCS expansion exhibit much more 
complicated behavior than the straightforward power counting in e would indicate. 

The study in KSMK1, however, is subject to elaboration in at least two points. 
The first is that, for the nucleon propagator, KSMK1 used the standard HB^PT 
propagator instead of the MCS propagator (see below), and the consequences of this 
eclectic treatment need to be investigated. The second point is that, in KSMK1, only 
the nucleons and pions appear as effective degrees of freedom, with the A-particle 
field integrated away. In view of the rather high momentum transfers involved in 
the NN — > NNn reactions, the influence of including the A-particle as an explicit 
degree of freedom is well worth examination. Since however this second point has 
been discussed in Refs. [8j [12] H we concentrate here on the first point. Furthermore, 
we limit ourselves here to the study of the transition operators (or equivalently the 
transition amplitudes in plane-wave approximation) even though the actual transition 
amplitudes need to be calculated with the use of the distorted-waves (DW) for the 
initial and final two-nucleon systems. Previous studies (see, e.g., [U El [8], [13]) have 
shown that, in order to reproduce the experimental data, it is crucially important to 
evaluate higher chiral-order contributions to the transition operator. In the present 
work, therefore, we focus on the convergence property of the MCS expansion of the 
transition amplitudes calculated in plane-wave approximation. As mentioned, the 
NN^NNtc reactions serve as a pilot case for nuclear reactions that involve rather 
high energy-momentum transfers, and it is hoped that our present study (despite 
its stated limitations) will shed some light on the general issue of an effective-field 
theoretical treatment of those reactions. 

Of the various possible isospin channels for the NN^-NNw reactions, we con- 
centrate (as we did in KSMK1) on the pp-^ppn reaction. The amplitudes for the 
iViV — > NNn reactions are in general dominated by a pion-rescattering diagram in 

2 In Ref. [8] it is convincingly argued that at next-to-leading-order (NLO) in MCS, the diagrams 
involving a A cancel for s-wave pion production, i.e. A will only enter at higher orders in MCS. 



2 



(a) Resc (b) VII 



Figure 1: One-pion-exchange rescattering diagrams of order e 3 (in MCS) for the 
threshold pp-+ppir° reaction. Diagram 1(a) represents the lowest-order non- vanishing 
one-pion-exchange contribution (NLO in the W-scheme) , while diagram 1(b) repre- 
sents its recoil correction (NNLO in the W-scheme). Adopting the same convention 
as in Ref. [15] for labeling the diagrams, we call diagram 1(a) the rescattering dia- 
gram, or the Resc diagram for short, and refer to diagram 1(b) as diagram VII. 3 The 
re-scattering vertex in the Resc diagram comes from £S = Ci N^Tt X+ )N + 
whereas diagram VII contains the re-scattering (recoil correction) vertex of the next 
chiral order lagrangian C^. 



which rescattering is caused by the Weinberg- Tomozawa (WT) interaction term, but 
this particular diagram does not contribute to the pp — * ppn reaction. This feature 
makes the pp-^pprc reaction uniquely suited for investigating the behavior of higher 
order terms in s-wave pion production^] 

We have derived in Ref. [15] all the possible two-body transition operators for 
s-wave pion production in the pp—yppn reactions up to next-to-next-to-leading order 
(NNLO) in the W-schemeQ The diagrams generating these operators can be catego- 
rized into several groups: the pion-rescattering diagram illustrated in Fig. 1(a) and 
its recoil correction diagram (diagram VII) shown in Fig. 1(b); the two-pion exchange 
(TPE) diagrams (diagrams I~IV) in Fig. 2, the vertex loop-correction diagrams (di- 
agrams V and VI) in Fig. 3, and the contact-interaction diagram (CT diagram) in 
Fig. 4@ In what follows, the diagram in Fig. 1(a), which specifically represents the 



3 It should be mentioned, however, that the s-wave pion production amplitude represents a more 
complicated case than the p-wave pion production amplitude; the convergence property of MCS for 
the latter has been discussed in Ref. [14] . 

4 The structure of the one-body transition operators is well known in HBxPT. 

5 In labeling the diagrams in Figs. 1~4, we are following the same convention as in Ref. [15]; this 



I 



II 



III 



IV 



Figure 2: The two-pion-exchange (TPE) diagrams considered to be of order e 2 in 
the formal counting in MCS (NNLO in W-scheme). We only show one representative 
diagram for each group, suppressing similar diagrams belonging to the same group 
(see Ref. [15] for details). The diagram labels, I, II, III and IV, follow the convention 
used in Ref. [15] . All the vertices arise from the lowest chiral order lagrangians, C^ N 
and C®>, see Ref. [TBI - 




Figure 3: One-pion-exchange diagrams with vertex correction belonging to order e 5 in 
MCS (NNLO in the W-scheme). The diagram labels, V and VI, follow the convention 
used in Ref. [T5] . 




CT 

Figure 4: The counter-term (CT) diagram of order e 3 (NNLO in the W-scheme). 

Table 1: The orders of diagrams in the W-scheme and MCS. The labels of the dia- 
grams conform with those used in Ref. [15]. When a diagram is assigned 0(e n ) in 
MCS, it means that its "leading term" is 0(e n ) with the understanding that there 
can also be terms higher order in e ( "sub- leading terms"). 



Diagram 


Resc 


I 


II 


III 


IV 


V 


VI 


VII 


CT 


W-scheme 


NLO 


NNLO 


NNLO 


NNLO 


NNLO 


NNLO 


NNLO 


NNLO 


NNLO 


MCS 


0(i 3 ) 


0(e 2 ) 


0(e 2 ) 


0(e 2 ) 


0(e 2 ) 


o(i 5 ) 


0(~e 5 ) 




o(i 3 ) 



lowest-order non- vanishing one-pion-exchange contribution to pp — > ppir , shall be 
simply referred to as the rescattering diagram, or the Resc diagram for short. Tabled] 
shows the power counting of these diagrams (Figs. 1 ~ 4) in the W-scheme (pow- 
ers in e) and in MCS (powers in e). The table indicates that the two schemes give 
significantly different classifications of these diagrams. Diagrams I ~ VII and the 
CT diagram, which are all categorized as NNLO terms in the W-scheme, belong to 
different orders in MCS, ranging from 0(e 2 ) to 0(e 5 ). Furthermore, the Resc term, 
which in the W-scheme is of NLO and thus of lower chiral order than any other terms 
in the table, belongs to (9(e 3 ) in MCS, and thus of higher order in e than Diagrams 
I ~ IV. These features invite us to investigate the actual numerical behavior of the 

leads to a somewhat awkward situation here that the diagrams are not numbered according to the 
order of their appearance in the text; "diagram VII" shows up earlier than diagrams I, II, III, etc. 
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diagrams listed in Table [T] in the context of MCS and examine whether MCS indeed 
provides a useful guide in organizing higher order contributions. The main purpose 
of this article is to report on such a numerical investigation. 

A major difference between the W-scheme [8] and MCS jH] is associated with the 
different treatments of the nucleon propagator (see e.g. Ref. [IS]), and we briefly 
recapitulate this point here. In HB%PT the nucleon momentum is written as 
Qt 1 = rriN v 11 + q fJ ', where rriN is the nucleon mass, and is the four- velocity, which 
may be chosen to be = (1; 0); it is assumed that |g^| <mjy. In terms of the 
Feynman propagator for the heavy nucleon is expressed as 

S N ( q ) = i $ + niN =i ( v . q + q2 \ 1 ( 1 + 70 _|_ 7o ffo _ 

Q 2 — m 2 N \ 2rriN ) \ 2 2m^ 2m^ J 

In one version of non-relativistic HB^PT, the free heavy-nucleon Lagrangian is chosen 
in such a manner that the heavy-nucleon propagator is given by 

S N (q) = — , (2) 
v ■ q 

and the difference between the propagators in eqs. (|IJ) and (j2|) is treated as perturbative 
recoil corrections which are accounted for in higher chiral order Lagrangians. For 
convenience, we refer to this approach as the W-schemeJ§ The work in Refs. [TH [15] 
is based on the use of Sx(q) given in eq.([2J). 

Meanwhile, in MCS, one rearranges the expression in eq.flTJ to adapt it to a specific 
kinematic situation pertaining to the NN^NNtt reaction at threshold. In order for a 
threshold pion to be produced from two nucleons, the four-momenta of the incoming 
on-shell nucleons 1 and 2 (in the CM system) with m N v^ subtracted must be pi = 
(m n /2,p) andp 2 = (m n /2, —p), where \p\ = ^m N m n +ml/4: = % /m Jv m 7r [l+e 2 /8+0(e 4 )]. 

(Recall e = \Jrn^Jm N .) In what follows, for notational simplicity, we write p instead of 
P\. When an incoming on-shell nucleon of four- momentum p produces a pion of four- 
momentum I, the off-shell nucleon propagator has a momentum p—l, where the loop- 
integral four-momentum / behaves like: v-l — Iq ~ Ip^^^/WiNTn^, \l\ ~ \p\ ^^/miym^. 
Let the nucleon propagator in eq.(0Q) be rewritten as 

s N{p -i)J v . (v - o - (g=Sl + I'-fe - Ol'Y V l±7._r(gzi) + T,fa - h)\ 

KF I 2m N 2m N M 2 2m N 2m N I 

6 It should be remarked, however, that Weinberg [17] discussed problems associated with using 
the static nucleon propagator, Eq. l[2|). in the two-nucleon systems, and that he proposed a possible 
remedy, which however upset the original counting scheme. 
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i —v-l+v-p- 



P ' 



P-21-p , \po-U' 



2m 



N 



2m 



N 



2m 



l+7o 1-(P-1) lo(Po-lo) 



N 



2rriN 



2mjy 



Since v-p = ^rf and v-l~ ^miym n , the difference 



v-p 



P ' 



2m 



N 



2 



777^777^(1 +e 2 /4) 



2m 



e 3 



(3) 



(4) 



appearing in eq.Q is of higher order in e compared to v ■ I and the other terms in 
eq.Q and therefore can be dropped. This treatment of the p 2 /2m N term is to be 
contrasted to that in the W-scheme, where the p 2 /2m^ term is included in the higher 
chiral order lagrangian and treated as a perturbative correction. We now reorganize 
eq.Q as 



S N (p-l)=i 



-vl- 



l 2 -2l-p (p -l 



2m 



N 



2m 



N 



l+7o 7-(p- , 7oG»o - ^o) 



2m N 



2m 



+ 



N 



(5) 



Here, compared with the leading term v-l, the terms in the curly brackets are of order 
e or higher. Expanding eq.© in powers of l/2mjv we arrive at the expression given 
in appendix B of Ref. [18] : 



s N (p-i) 



vl 



l+7o 



l 2 -2l-p 
2m N l 



7*(P- , ( v ' 1 ) A-7o 



2m 



N 



2m 



N 



+ 



(6) 



where the three terms containing the factor l/2mjv are of order e relative to 1. In 
MCS, one applies the lowest-order non-relativistic approximation to eq.Q (letting 
7o — > 1 and ignoring the e corrections) and adopts as the heavy-nucleon propagator 



S N (p-l) 



vl 



(7) 



In Ref. [18] . the \j2m^ term within the curly brackets in eq. ([6]) was classified as a 
heavy-nucleon propagator recoil correction, while the other two \j2m^ terms were 
classified as higher order vertices in the the Lagrangian. In the present work, following 
Ref. [18J, we use eq.© without considering these corrections. 

We investigate numerical differences in the transition amplitudes resulting from 
the use of the two different expressions, eqs.© and (j7]), for the nucleon propagator. 



7 



As mentioned, the analytic expressions for the pp^ppn transition operators were 
calculated in Ref. [T5] up to NNLO in the W-scheme [i.e., the HB^PT propagator, 
eq. (jSJ) was used]0 It turns out that the amplitude expressions in MCS can be readily 
obtained from those given in Ref. [15]. The dependence on v-p in the amplitudes 
enters only through the HB^PT propagator eq. (J2]) with q = p — I. This implies that 
the amplitudes resulting from the use of the MCS nucleon propagator, eq. (J7J, can be 
simply obtained by formally replacing v-p with zero (v-p—*0) in the expressions given 
in Ref. [15]. These expressions are valid for arbitrary kinematics but, for the present 
purpose, we simplify them with the use of the fixed (frozen) kinematics approxima- 
tion (FKA), wherein the external energies and momenta appearing in the particle 
propagators and vertices are "frozen" at their threshold values. 

In FKA, the operator corresponding to each of the diagrams in Figs. 1~ 4 are 
written as 




where p (p r ) is the relative three-momentum in the initial (final) pp state {p\—p2 = 2p, 
p{ —P2 = 2p'), k = p — p', x = p-p', and S = — 02). The partial- wave projected 
form of t(p,p',x) is written as [13J 

J W = - ( Sn ^ J J P 2 dp P ' 2 dp'J dxi[)i So (p') t(p,p',x)( y p-p'x)iJip (p) (9) 

where ip a (p) is the a partial- wave ( 1 iSo for the initial state and 3 Pq for the final state) 
For the t(p, p', x)'s corresponding to the two-pion exchange (TPE) diagrams (diagrams 
I, II, III and IV in Fig. 2), it is informative to decompose each of them into terms 
with definite asymptotic /^-dependence as [HI [19] 

t(p,p',x)= t 1 (g A /(8f n )) 2 \k\+t 2 (\n{\k\ 2 /A 2 }) + t 3 + 5t(p,p',x), (10) 

where, in the limit of k —>■ 00, £3 is k- independent, and 5t(p,p',x) is (9(/c _1 ). The 
analytic expressions for tj's (i — 1, 2, 3) can be extracted [19] from the amplitudes T 
given in Ref. [15]. The first term with t\ in eq. ffTOl) is the "leading part" of 0(e 2 ) in 
MCS, which was evaluated in Ref. [7]; the remaining "sub- leading" terms in eq. ffTOj) 

7 The amplitude for diagram II given in Ref. [15] has the wrong sign. 

8 In this work we adopt the plane- wave approximation, which corresponds to the use of the wave 
functions of the generic form ip(p) = 5{p — p n)/p 2 , where p on is the asymptotic on-shell nucleon 
momentum. 
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were not considered in Ref. [7J. For the reason to be explained later, we introduce 
the "sub- leading" terms of 0(i 3 ) as 

t*(p,p',x) = t(p,p',x)-t 1 (g A /(8f^) 2 \k\ (11) 
and, correspondingly, 

= _ (^^) l o °°p 2 dpp' 2 dp' J'dx ^ So (p') t*(p,p',x)(p-p'x)^ Po (p) (12) 

We remark that, for diagram I, we have t*(p,p',x) = t(p,p',x) and J[t*] = J[t], 
since the t\ term exists only for diagrams II, III and IV [Til [19] . In what follows, we 
investigate the numerical behavior of J[t] and J[t*]. For the sake of definiteness, we 
concentrate here on a representative near-threshold case where the kinetic energy of 
the incident proton (in the laboratory system) is 7} a & = 281 MeV. 

We first consider the two-pion exchange (TPE) contributions, coming from dia- 
grams I ~ IV depicted in Fig. 2. Since the leading parts of the TPE diagrams of 
0(e 2 ) are known to cancel among themselves [TJ [11] Jj our main concern here is the 
behavior of t*(p,p',x), eq. ljTTI) . or equivalently J[t*], eq. (fT2l) . We show in Table[2]the 
values of J[t*] corresponding to each of diagrams I ~ IV. In the last column we also 
give J[t]Resc, which is the contribution to J[t] from the rescattering diagram (Resc) 
(Fig.l). The second row gives J[t*] calculated in MCS, while the third row shows 
J[t*} calculated in the W-scheme [llj. Since the rescattering diagram contains no 
nucleon propagators, the W-scheme and MCS give the same value for J[t]R esc • The 
fourth row gives J[ti] representing the contribution of the leading t\ term; we remark 
that, for J[ti], both MCS and the W-scheme give the same result. For comparison, 
we present in the fifth row the value of J[t] calculated in Ref. [TH] with the use of the 
W-scheme. 

Comparison of the J[t*; MCS\s and J[t*; W]'s in Table [2] indicates that the dif- 
ferent treatments of the heavy-nucleon propagator between the W-scheme and MCS 
affect rather appreciably the individual values of J[i*]'s for diagrams I ~ IV. Thus, 
for a formally consistent check of MCS it is important to use the MCS nucleon prop- 
agator, and in what follows we shall mostly concentrate on the results obtained with 
the MCS nucleon propagators. 

According to MCS, J[t*] for the TPE diagrams should be of 0(e 3 ), and J[t]R esc 
should also be of C(e 3 ), see Table [2J Thus, if MCS is a reasonable counting scheme, 

9 This can also be seen in our Table [2j Since t\ = : — 1 : — 1/2 : 3/2 for diagram I, II, III and 
IV, respectively, the J[ti]'s in the table are found to be in the ratio of : —2 : — 1 : 3 and they add 
up to zero. 
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Table 2: The pp^ppn amplitudes J, eqs. ([9]) and ( fT2j) . calculated for Ti a b = 281 MeV 
in the plane- wave approximation, and in the frozen kinematics approximation (FKA). 
The labels, I~IV, correspond to the diagrams I~IV depicted in Fig. 2, while the col- 
umn labeled "Sum" gives their combined contributions. The last column shows the 
contribution to J[t] from the lowest-order non- vanishing one-pion-exchange rescat- 
tering diagram, the Resc diagram (Fig. 1(a)). The second row shows J[t*] calculated 
in MCS and formally of 0(e 3 ), while the third row gives J[t*] evaluated in the W- 
scheme [15J . The fourth row gives J coming from the leading t\ term of MCS order e 2 ; 
since J[ti;MCS] = J\t\\ W] we simply write J[ti]. The fifth row gives J[t] calculated 
in Ref. [T5] with the use of the W-scheme. J[t* W] + J[h] = J[t; W], see eq.lfTTl). 



Type of diagram 


I 


II 


III 


IV 


Sum 


Resc 


J[t* MCS] x 10 2 


-1.4 


-13.2 


-23.1 


12.5 


-25.2 


8 


J[t* W] x 10 2 


-5.6 


-6.5 


-29.9 


4.9 


-37.1 


8 


J[h] x 10 2 





-46.7 


-23.3 


70.0 







J[t; W] x 10 2 


-5.6 


-53.2 


-53.2 


74.9 


-37.1 


8 



J[t*] and J [t] R esc are expected to be of comparable magnitudes. A comment is in order 
here, however, on the meaning of "comparable magnitudes" . Usually, two quantities 
that differ by a factor of 2~3 are considered to be of comparable magnitudes (or of 
the same order) but, in the present case where the expansion parameter is e ~ 1/3, a 
difference by a factor of ~3 may be interpreted as representing a different order in e. 
Ideally speaking, one could make the situation clearer by comparing terms that differ 
by two orders in e. In our present study, however, this is possible only for certain 
classes of diagrams (these cases will be discussed later in the text). For the other 
cases in which MCS can be tested only within one order in e, we shall adopt the 
following viewpoint: If a deviation from the behavior expected from MCS lies within 
a factor of ~3, we categorize it as reasonably consistent with MCS, in the sense that 
the deviation does not constitute definitive evidence for a breakdown of MCS. 

In looking at Table we first concentrate on diagrams II, III and IV, leaving out 
diagram I for a while. We note that J[t*; MCS] for diagrams II and IV are close to 
-^Resc, while J[t*; MCS] for diagram III is larger than J[t]R esc by a factor of ~30 
Thus the numerical behavior of the J[t*\ MCSj's for diagrams II, III and IV is either 

10 In what follows, we are mostly concerned with the absolute value of J rather than J itself but, 
for the sake of simplicity, we shall refer to the absolute value of J simply as J (when there is no 
danger of confusion). 
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Table 3: The pp^ppir amplitudes J[t], eq. (jHJ), calculated for Tz Q & = 281 MeV in 
the plane-wave approximation and in FKA. The labels, V, VI and VII, correspond to 
the diagrams V, VI and VII depicted in Figs. 3 and 1, respectively. The last column 
gives the contribution to J[t] from of the rescattering diagram (Fig.l). The second 
and third rows show, respectively, J[t] calculated in MCS and in the W-scheme [TB] . 



Type of diagram 


V 


VI 


VII 


Resc 


J[t; MCS]xl0 2 


3.4 


-2.3 


20.8 


8 


J[t; W] x 10 2 


1.4 


1.1 


20.8 


8 



consistent or reasonably consistent with MCS. We find a similar situation in compar- 
ing J[t*; MCS] (second row) and J[ti] (fourth row). According to MCS, J[t*] should 
be one order higher in e than J[ti\. For diagrams II and IV, J[t*; MCS]'s exhibit a 
slight over-suppression (beyond e ~l/3) but they are still reasonably consistent with 
MCS. For diagram III, J[t*; MCS] ~ J[ti] and hence no suppression is seen, but the 
level of deviation from the MCS rule is again within a factor of ~3, presenting another 
reasonably consistent case. 

According to Table [2} J[t*; MCS] for diagram I is smaller than J[t]R esc by a 
factor of ~6, in contrast to the other TPE diagrams. It is to be noted however that, 
for diagram I, the leading term itself exhibits an unusual behavior; viz., although a 
straight-forward MCS counting indicates that diagram I, containing only one nucleon 
propagator, should be of order e 2 , the actual calculation shows t\ = [TT]J"1 It is 
thus likely that some (possibly accidental) extra suppression mechanism is at work for 
diagram I. Since the existence of this type of extra suppression does not necessarily 
signal a breakdown of a perturbation series, we may take the view that the results 
for diagram I is essentially consistent with MCS. Based on the above discussion, 
we conclude that the overall behavior of our numerical results shown in Table [2] 
is reasonably consistent with MCS. As mentioned, we are ignoring here the recoil 
correction to the MCS propagator eq. (j7j). We will return to this issue later in the 
text. 

We next discuss the behavior of the other higher-order contributions coming from 
diagrams V and VI (Fig. 3) and diagram VII (Figs. 1(b)). In diagram VII the pion- 
nucleon rescattering vertex is given by C^, see eq.(C.3) in Ref. [20]; this vertex 
contains the recoil correction to the pion-nucleon vertex in the lower chiral order re- 

11 This corresponds to the remark made in Ref. [7] that diagram I is "beyond next-to- leading 
order" . 
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scattering amplitude (Resc) in Fig. 1(a). In Table [3] we show J[t] for diagrams V, VI 
and VII, calculated in MCS (second row) and in the W-scheme (third row); this latter 
has been taken from Ref. [15]. According to MCS [Sj, diagrams V and VI, which are 
the pion s-wave re-scattering diagrams containing a pion loop at one nucleon vertex, 
should be suppressed by e 2 ~ 1/7 compared to diagram VII. The numerical results 
in Table [3] are in conformity with this expected suppression. Meanwhile, diagrams 
VII and Resc (shown in Fig.l) should both be of order (e) 3 in MCS 0, [8]. Table [3J 
indicates that J[t] for diagram VII is larger than J[t]R eS c by a factor of ~3, but here 
again this level of deviation from the MCS prediction is regarded as constituting a 
reasonably consistent case. It is informative to study in more detail the origin of the 
differences between these two contributions. For threshold kinematics (FKA), the 
pion-nucleon re-scattering vertices in the Resc diagram and diagram VII are given, 



Formally, these two vertices do seem to be of the same order in MCS. However, the 
well-known smallness of the irN iso-scalar scattering length requires a substantial 
cancellation between the three q terms in eq. ffTBl . This particular situation causes 
the expression in eq. (jl~4]) to become larger than that in eq. (|T3|) by almost a factor of 3. 
These features illustrate intrinsic subtlety we encounter in discussing the convergence 
property of MCS beyond the level of a factor of ~3. On the other hand, according 
to the W-scheme, the Resc diagram belongs to NLO and diagram VII to NNLO (see 
Table 1). Our numerical results exhibit a definite deviation from the prediction based 
on the W-scheme. 

We have seen above that, despite the fact that the expansion parameter e in MCS 
is not very small (e ~ 1/3), MCS is likely to provide a useful guide in organizing 
the higher-order terms. In this connection, it seems worthwhile to examine J[t] as 
a function of e by allowing e to vary from its physical value to smaller values. This 
may offer information somewhat different from that obtainable by comparing the 
magnitudes of the individual diagrams at the fixed value of e. We explore here how 
the contributions of the diagrams behave when the pion mass m n is artificially varied. 
To extract the non-trivial m^-dependence of J, we define the plane-wave amplitude 

12 For the vertex for the Resc diagram, see eq.(A.29) in Ref. [16]; for the vertex for diagram VII, 
see eq.(C.3) in Ref. [20]. 



respectively, byr^l 



Resc 



VII 




(13) 



(14) 
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X by 



m = -(^) 2|*| A' (15) 

As for the "leading term" of the TPE amplitude (viz., the ti-term contribution), we 
take as a representative the J[ti] for diagram II and introduce 

J[h; diagram II] = - 2\k\ X tl (16) 

which gives (ti = 1 for diagram II) 

X tl = (g A /(8f 2 n )) 2 \k\. (17) 

In Table H] we show as functions of the values of X for individual diagrams 
calculated in MCS. The graphical representation of X/m n as a function of is 
given in Fig. 5. We expect from eq. ffTTl) that X tl ~ y/m n for — > 0, and the 
numbers in the second column of Table H] indicate that this expectation is borne out. 
According to MCS, the A's for diagrams II, VII and the Resc diagram are expected 
to be linear in m n for m n — > 0. The numerical results in Table H] are in agreement 
with this MCS expectation. According to Fig. 5, the As for the TPE diagrams I, 
III and IV appear to have a more complex m n dependence than implied by a naive 
application of MCS. The As here diminish less rapidly for m n — >0 than the expected 
linear dependence on m n . We remark that in MCS the chiral log factor, h^m,,-), 
which arises from the loops of the TPE diagrams, is treated as a constant factor that 
does not affect the MCS power counting. Meanwhile, the existence of the t2-term in 
eq. ffTUl) implies that the m^-dependence in the chiral log can disturb the simple power 
counting in MCS, and Fig. 5 indicates the chiral log behavior is dominant at low 
values. Finally, according to MCS, the As for diagrams V and VI are expected to 
vary as m 2 for — > 0. Our numerical results indicate that they go to zero faster 
than m n but somewhat slower than m\. 

Besides the diagrams so far considered, the five-point vertex counter-term diagram 
(CT diagram) shown in Fig. 4 is also of order e 3 in MCS [S] , and we now give a brief 
discussion of the numerical behavior of the CT diagram. An aspect that distinguishes 
this diagram from the others is that the LEC here is an unknown parameter whereas 
the LECs in the other diagrams are predetermined from separate sources of infor- 
mation. In the absence of experimental data needed to determine the relevant LEC, 
we rely here on the resonance saturation prescription and use the a- and u;-exchange 
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Table 4: The plane-wave amplitude X [eg. (11 51) ] (in units of 10~ 2 fm 3 ) calculated for 
Ti ab = 281 MeV in FKA; e.g., the second column gives the values of X tl , eq. (|T7|) . For 
the cases where MCS and the W-scheme give different results, MCS is used. 



Diagram 


ti 


£*: I 


t*: II 


£*: III 


£*: IV 


t: Resc 


t: V 


t: VI 


t: VII 


m« = 140 MeV 


92.6 


2.7 


26.2 


45.9 


-24.8 


-15.7 


-6.8 


4.6 


-41.5 


100 


77.9 


1.5 


18.6 


43.8 


-19.5 


-11.6 


-3.8 


2.6 


-29.9 


50 


54.7 


0.3 


9.4 


32.6 


-11.3 


-6.1 


-1.1 


0.8 


-15.1 


10 


24.3 


-0.1 


2.1 


10.9 


-2.9 


-1.3 


-0.06 


0.04 


-3.1 



0.05 

0.04 

5- 0.03 
E 

h — _ 

E* 0.02 
X 

0.01 






1x10 


i 


II x 10 

III 


\ 


IVx(-10) 










— ^ . 




\ 







20 40 60 80 100 120 140 



[MeV] 

Figure 5: For each of the TPE diagrams, J/m T is plotted as a function of m^; the 
amplitudes X are defined in eq. (|T5l) . 
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mechanism considered in Lee and Riska's work [3]. In this model, the amplitude T 
[eq.flS])] for the CT diagram is given by, see e.g. Ref. [5] 



T, 



CT 



9A^g 



f(4- 


h4) 




mi) 



•2(7 _|_ JUJ j ^ ^) 



mi 



[1 + C U ) {ff 1 xa 2 )-k 



where u q = \jm\ + q 2 is the energy of the outgoing pion and P = p + p 'o A 
caveat here is that, since the "a" exchange is considered to represent a scalar part of 
correlated two-pion-exchange in the AW potential, and since we have already taken 
into account some TPE diagrams, there is a danger of double-counting if we include 
the entire "a" exchange contribution to the CT diagram. Without going into this 
issue, we present here the individual contributions of the "<r"- and cu-exchange and 
compare them with the contributions of the other diagrams that are expected to be 
of order e 3 in MCS. If again we restrict ourselves to threshold kinematics (p ' = and 
k — P — p), eq.l fTSl can be effectively rewritten as 



T, 



CT 




4-+4-[i+2(i+c w )] 



mi 



mi- 



(19) 



which, with the use of £/(4>ir) = 5, gl/(An) = 10, 7^^ = 600 MeV and m w = 780 MeV, 
leads to 



T CT = i^A S-p{l3.8GeV- 3 + 16.4GeV- 3 [l + 2(l + CJ]} 



(20) 



Referring to the relation between T and J, see eqs. (jHJ) and (jHJ), we estimate Jqt from 
eq.(|20l) to be 



J CT = -0.052 - 0.062 [3 + 2C„ 



(21) 



where the first (second) term comes from u-exchange (tu-exchange) The magnitude 
of the Jqt amplitude is to be compared to J(t*] MCS) in the second row in Table El 
Within the factor of ~ 3 the a- and tu-exchange contributions are comparable to the 
C(e 3 ) TPE amplitudes in Table 2. Although, as stated earlier, there is a question of 



13 Between the spin singlet and triplet states (between 1 5'o and 3 Po), the operator i (a\ x a?) is 
equivalent to (— 2)S. 

14 In the literature we find \CJ\ < 1 but its value is not well determined; see, e.g., Refs. [2T| [22] . 
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possible double-counting concerning the "a" -exchange contribution, it is to be noted 
that the CT amplitude coming from tu-exchange alone is of expected MCS magnitude. 

To summarize, we have calculated contributions to the threshold pp — > ppn re- 
action amplitude with the use of the MCS heavy-baryon propagator, studied the 
convergence property of the expansion based on MCS, and compared the results with 
those obtained in the W-scheme. We have considered all the diagrams up to NNLO in 
the W-scheme and, taking advantage of the fact that these diagrams belong to widely 
different orders in MCS (ranging 0(e 2 ) to 0(e 5 )), we have examined whether the 
actual numerical behavior of these diagrams substantiates the MCS predictions. Our 
numerical results indicate: (a) The terms that are expected to be of the same order 
in MCS follow that pattern within a factor of ~ 3; (b) The terms that are expected 
to differ by one order in e in MCS also follow that pattern within a factor of ~ 3; (c) 
The terms that are predicted by MCS to be suppressed by two orders in e (e 2 ~ 1/7) 
clearly show the expected suppression. Thus all the diagrams studied here exhibit 
numerical behaviors that are (in our definition) reasonably consistent [Cases (a) and 
(b) above], or consistent [Case (c)]. These results lead us to conclude the following 
(i) Even though the expansion parameter e in MCS is not very small (e ~ 1/3), MCS 
provides a useful semi-quantitative guide in organizing the higher-order terms; (ii) 
MCS is superior to the W-scheme in organizing the relative importance of the higher 
order terms. This last conclusion supplements the finding in Ref. [2], which shows 
that MCS converges for p-wave pion production. 

A better test of the convergence property of MCS would be to calculate next-order 
contributions for all the diagram under consideration and compare the terms that are 
expected to differ by two orders of e, but this is beyond the scope of the present work. 
As mentioned we have neglected the 0(e 3 ) contributions to the reaction amplitude 
arising from the recoil correction to the MCS propagator, eq. ((?]), and we have also 
ignored the additional vertices discussed in connection eq.flTj). We remark, however, 
that these contributions do not directly affect the behavior of the contributions of the 
diagrams studied in the present investigation. 
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